Intermediate quantum maps for quantum computation 



in 
o 
o 

(N 

Oh- 
< 



> 
O 
m 

(N 

O 
in 
o 

^— > ■ 

G ■ 
cd 

^ : 
cr. 



0. Giraud and B. Georgeot 
Laboratoire de Physique Theorique, UMR 5152 du CNRS, 
Universite Paul Sabatier, 31062 Toulouse Cedex 4, France 
(Dated: April 29, 2005) 

We study quantum maps displaying spectral statistics intermediate between Poisson and Wigner- 
Dyson. It is shown that they can be simulated on a quantum computer with a small number of gates, 
and efficiently yield information about fidelity decay or spectral statistics. We study their matrix 
elements and entanglement production, and show that they converge with time to distributions which 
differ from random matrix predictions. A randomized version of these maps can be implemented 
even more economically, and yields pseudorandom operators with original properties, enabling for 
example to produce fractal random vectors. These algorithms are within reach of present-day 
quantum computers. 
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The study of quantum information has attracted more 
and more interest from the scientific community in the re- 
cent past. Quantum communication and quantum com- 
putation have been shown to be deeply different from 
their classical counterparts. Algorithms have been built 
showing that quantum computers can outperform classi- 
cal devices for some problems Q. In particular, quantum 
mechanical systems can be simulated much faster on a 
quantum computer An especially congenial class of 
models corresponds to quantum maps. They can dis- 
play complex dynamics, but can be described by simple 
evolution operators. It has been shown that the quan- 
tum baker's map Q, the quantum kicked rotator U, the 
quantum sawtooth map [5j can be simulated efficiently 
on a quantum computer. Still, even the simplest mod- 
els are difficult to implement on the small-size quantum 
computers experimentally available, and only the quan- 
tum baker's map has been implemented to date with 
three qubits |(J. General algorithms have been proposed 
to probe phase space distributions Q, fidelity decay ||, 
form factors or localization length 01 01 such maps. 
Quantum maps have been also used as testgrounds to 
study the production of entanglement in quantum sys- 
tems ^lj- Very recently, they have been used as models 
to build pseudorandom operators which can be efficiently 
implemented on quantum computers [f^ ]. 

The models envisioned so far correspond mostly to 
chaotic maps, where the most complex behavior is ex- 
pected to manifest itself, and eigenvalue statistics are 
close to those of the Wigner-Dyson ensembles of Random 
Matrix Theory (RMT). Some integrable maps, which are 
expected to follow Poissonian statistics, have been also 
studied. However, it has been observed recently that 
some systems obey intermediate statistics, characterized 
by a level repulsion and a Poisson-like behavior at long 
distance in energy spacings. They can be modeled by the 
so-called semi- Poisson statistics for which correlation 
functions can be explicitly calculated, giving energy level 
spacing distributions Pp{s) cx s /3 e~^ +1 - )s . Such distri- 



butions have been first observed at the Anderson metal- 
insulator transition for electrons in disordered systems 
[l3| , and later in pseudointegrable systems 0] or in cer- 
tain diffractive billiards 0] . They are usually associated 
with fractal properties of eigenstates. 

In this paper, we study from the viewpoint of quantum 
computation a one-parameter family of quantum maps 
recently introduced | l6| . They can be expressed in a par- 
ticularly simple way, yet the spectral statistics display a 
wide range of semi-Poisson distributions depending on 
the value of the parameter. We first show that such 
maps can be simulated efficiently on a quantum com- 
puter, and can be implemented with a remarkably small 
number of qubits and quantum gates. In particular, it 
represents an ideal testground for algorithms proposed 
in recent years which aim at measuring the form factor 
or the fidelity decay using only one qubit of quantum 
information. The properties of these maps being inter- 
mediate between chaos and integrability, we study how 
this translates in the matrix element distributions and 
the entangling power. At last, we show that a certain 
generalization of these maps which corresponds to a new 
ensemble of random matrices with semi-Poisson statistics 
recently proposed in |l7| can be implemented even more 
efficiently on a quantum computer. They can be used 
as a way to produce a family of pseudorandom operators 
with new properties related to fractal behavior. We note 
that different distributions interpolating between Poisson 
and Wigner-Dyson, built from a partial randomization 
procedure, were studied in [lij . 

We start with the classical map defined by p = p + 
a (mod 1) ; q = q + 2p (mod 1) where (p, q) is the 
pair of conjugated momentum (action) and angle vari- 
ables, and the bars denote the resulting variables af- 
ter one iteration of the map. The quantization of this 
map yields a unitary evolution operator which can be 
expressed in momentum space by the N x N matrix 

TT , - c*P(-2»7rp' 2 /jV)(l-cxp(2»7rjVa)) [7TI alfprnflti vel v 

u pp' ~ i-cxp(2iK(p-p>+Na)/N) l±9> or alternatively 
in operator notation U = e~ 2l7T P ' N e 2l7Ta ^ . For generic ir- 
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FIG. 1: Nearest-neighbor distribution of eigenvalues of U for 
N = 2 12 and from left to right a = 1/3, a = 1/5, a = (1 + 
v5)/2. Data are taken from the dissymmetrized spectra of 
a and 1 — a. Distribution corresponding to a — 1/5 (resp. 
a = (1 + \/5)/2) is shifted by 1 (resp. 2) along the s-axis. 
The solid curves from left to right correspond to semi-Poisson 
with /3 = 1/2, /3 = 3/2 and COE. 



rational a, the spectral statistics of U are expected to fol- 
low RMT (in this case COE). For rational a = a/b, a va- 
riety of different behaviors are observed |l6|. In particu- 
lar, it is conjectured from that for aN = ±l (mod b) 
the statistics of eigenvalues is of the semi-Poisson type 
Pp(s) with parameter (3 = 6/2 — 1. This can be checked 
in FigQ where the statistics of eigenvalues is plotted for 
different values of a (for N = (mod 4) there is an 
additional symmetry US — SU, where S qq > — (—l) q 5 qq >; 
each half of the spectrum should be considered separately 
jlfij). Thus U gives a set of quantum maps with statis- 
tics corresponding to natural intermediate distributions 
between Poisson and RMT, in a controllable manner. 

The map U can be implemented efficiently on a quan- 
tum computer. Indeed, the use of the Quantum Fourier 
Transform (QFT) allows to switch from position to mo- 
mentum representation. More precisely, for one iteration 
of U on a TV-dimensional Hilbert space with N = 2 nq , 
one first implements exp(2iTraq), which is diagonal in 
the ^-representation. This can be done efficiently using 
the binary decomposition of q: if q = J2j=o > then 
exp(2inaq) corresponds to the application of the n q one- 
qubit gates |0) -> |0), |1) -> exp(2«7ra2>')|l). Then by 
using a QFT one can shift from q to p representation, 
using n q {n q + l)/2 gates. In this representation, the sec- 
ond operator e - 2l ^P 2 / N j s diagonal. If p = ^j=o Pj^ '"> 
then exp(-2iirp 2 /N) = U 3u j 2 e^p{-2iiTp n p j2 2^ + ^ /N). 
To simulate it, one needs n 2 two-qubit gates applied to 
each qubit pair (ji, j 2 ), keeping the states |00), |01), 1 10} 
unchanged while |11> -> exp(-2in2^ + ^ /N)\ll) . Then 
a QFT brings back the wavefunction to the q representa- 



tion. In total, the evolution requires 2n 2 + 2n q gates to 
be implemented, of which 2n 2 q — n q are two-qubit gates. 
This is less than any other map proposed to date (includ- 
ing the sawtooth map except the quantum baker's 
map, which has already been implemented |(J. 

The algorithm above can be used as subroutine of other 
algorithms which aim at measuring quantum-mechanical 
quantities, with the attractive feature that it makes them 
very economical since the map evolution needs remark- 
ably little quantum gates. One can for example probe 
Wigner and Husimi phase space distribution functions 
, or investigate fidelity decay in presence of perturba- 
tion the fact that these maps correspond to spec- 
tral statistics intermediate between Wigner-Dyson and 
Poisson should translate into specific properties for these 
quantities. It has also been proposed to use a quantum 
computer to differentiate between quantum chaos and in- 
tegrability by evaluating the form factor at short times 
U In our case, the same algorithm can give much more 
information. The method in 0] adds one probe qubit to 
the system, performs U n , and uses two additional one- 
qubit gates to transfer the trace of U n to the probe qubit. 
Using such scattering circuits, real and imaginary parts of 
TiU n /N correspond to expectation values of Pauli opera- 
tors for the probe qubit. For COE or CUE, this quantity 
is of order 1/N for small n. In the case of intermediate 
statistics (a = a/6), one expects < Till 71 >„~ kVN, 
where the average is taken over the first iterates of U. 
The form factor at short times is then given by |k| 2 . 
In order to get the value of K with enough precision, 
one needs a number of quantum measurements of order 
iV (the number of values of n to average over depends 
only on 6 |lq | and does not vary with N). This im- 
plies a quadratic gain over classical computation. For 
N = (mod 4) , in order to dissymmetrize the spec- 
trum, one additionally has to perform the evolution of 
SU n (this only requires one extra controlled-phase gate), 
and the difference between the traces of U n and SU n 
gives the required quantity. Using this algorithm for large 
N enables to probe the form factor at increasingly short 
times and check the semiclassical conjectures 16J. The 
value of \k\ 2 also yields the level compressibility, which 
controls the behavior of the spectral number variance at 
infinity. It is possible that an exponential gain can be 
obtained through investigation of the fidelity decay 0. 
We note that the fidelity decay can be explicitly related 
to the form factor 19]. Both quantities can be probed 
using deterministic quantum computation with one sin- 
gle pseudopure bit 20], which together with the small 
number of gates needed should make these simulations 
very attractive for NMR quantum computation. 

It is instructive to study the entangling power of these 
maps. Indeed, entanglement is a key resource for quan- 
tum information 0. Quantum chaotic evolutions have 
been shown to generate entanglement distributions sim- 
ilar to the predictions of RMT ^lj- The entangle- 
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FIG. 2: (color online) Distributions of matrix element (top) 
and Q of column vectors (bottom) for iterates U n of U for 
N = 2 9 , in p representation. Matrix elements x are rescaled 
by y — N\x\ 2 . Only one symmetry class of S is taken. From 
lowest to topmost curve: black is a — 1/3, red is a — 1/5, 
green is a — (1 + a/5)/2, purple is CUE. Full lines are averaged 
over 10 5 < n < 10 5 + 1000, dashed line are averaged over 
1000 < n < 2000 (top) and 2000 < n < 3000 (bottom). 
Logarithms are decimal. 

ment can be measured by the average bipartite entan- 
glement between one qubit and the rest of the system 
Q = 2 — 2/n q Y^,k=i ^Pki where pk is the density oper- 
ator corresponding to the fc-th qubit after having traced 
out the rest |2l|. Fig [21 displays the distribution of matrix 
elements and Q for the column vectors of U n for large n. 
It is compared with the RMT prediction computed nu- 
merically from the parametrization method of [2^ . For 
irrational a, the spectral statistics follow COE, but both 
matrix element and Q distributions converge to CUE pre- 
dictions. For values of a with intermediate statistics, 
there is a convergence to a distribution different from 
CUE predictions. This holds in p representation: in q 
representation both quantities converge to CUE predic- 
tions (data not shown). In all cases the convergence is 
faster for matrix elements than for Q; the convergence 
rate is similar to what is obtained for chaotic maps \Ti\ . 

In |l2^ . pseudorandom operators inspired by quantum 
chaotic maps were built which efficiently create entan- 
glement and matrix element distributions close to RMT. 
In our case, it is possible to construct ensembles of 
Intermediate Statistics Random Matrices (ISRM) based 
on the quantum map U ^?|]- The procedure consists 
in replacing er 2 ™v / N in U by the diagonal matrix 
(in p representation) (e J * p <5 P p'), where <J> P are random 
variables cither independent ("non-symmetric case") 
or verifying $jv-p = &p ("symmetric case"). General 
arguments |l7| indicate that for a irrational eigenvalue 
statistics for ISRM follow COE (symmetric case) or 
CUE (non-symmetric case). For rational a = a/b, 



FIG. 3: (color online) Same as FigEJfor iterates of ISRM, for 
N — 2 s (non-symmetric case). Averages are made over 1000 
disorder realizations at n = 10 J (full lines), n = 1000 (dashed 
lines, top) and n — 2000 (dashed lines, bottom). 



and matrix size N — > oo obeying aN = ±1 (mod b) , 
eigenvalue statistics follow the semi-Poisson prediction 
with (9 = 6—1 (non-symmetric case) or/3 = 6/2—1 
(symmetric case). These ISRM can be implemented 
on a quantum computer, actually more economically 
than U. The only difference with the algorithm sim- 
ulating U consists in the replacement of e~ 2m v l N by 
(e^fSppi). This operator multiplies each basis state 
\p) by a Gaussian random phase. It can be simulated 
by choosing n q + n s independent and uniformly dis- 
tributed random angles (f>k, with 1 < k < n q , and <j)' k , 
1 < k < n s , for some integer n s . Applying the operator 

Ul=n s CNOT ifcJfe rifel! (^ fc (^)CNOT lfcJfc ) UlU MM 

multiplies each basis state by a random variable 
±<fii ± <f>2 ± • • • ± 4>' n -i±4>' n , which for large n s tends to 
a Gaussian random variable. Here Rj{4>) — exp^crj^) 
is the rotation on the j-th qubit by an angle 4>/2, 
CNOTjj (controlled-not) is the bit-flip on the j-th 
qubit conditioned by the i-th qubit, and the ik and 
jk are chosen randomly between and n q — 1. This 
transformation requires (3n s + n q ) gates instead of the 
n q gates needed for e ~ 2l7T P 2 / N . I n practice, n s is taken 
proportional to n q , and the simulation requires in total 
n 2 — n q + 2n s two-qubit gates and An q + n s one-qubit 
gates. This is quite close to what is needed to simulate 
the quantum baker's map already implemented [g, thus 
ISRM should be implementable as well in present-day 
quantum computers. 

The latter algorithm can be used to generate pseudo- 
random operators having original properties. Fig|3]shows 
that for ISRM matrix element and Q distributions con- 
verge to a limiting distribution which depends on a. For 
irrational a, the distribution corresponds to CUE predic- 
tions (even in the symmetric case where eigenvalue statis- 
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FIG. 4: IPR as a function of N = 2 n " for column vectors 
(circles) and eigenvectors (triangles) of ISRM (non-symmetric 
case), for a = 1/3, 1/5 (top), a = 1/7, (1 + V5)/2 (bottom). 
Straight lines show £ oc iV ' 80 , iV ' 93 , iV ' 97 , iV 100 (column 
vectors), and £ oc iV ' 54 , iV ' 69 , AT ' 76 , A^ 1 ' 00 (eigenvectors). 
Logarithms are decimal. 

tics follow COE, data not shown). For rational a where 
eigenvalue statistics follow semi-Poisson predictions, the 
limiting distributions are different from both COE and 
CUE distributions in p representation (in q representa- 
tion they converge to CUE predictions, data not shown). 
Intermediate eigenvalue statistics are also usually asso- 
ciated with fractal properties of eigenvectors. This was 
observed in the Anderson model at the metal-insulator 
transition or in pseudointegrable systems. As was seen 
in [J^ j this is also the case for the ISRM. In a quan- 
tum information setting, randomly chosen eigenvectors 
of ISRM can be obtained by using the phase estima- 
tion algorithm and measuring an eigenvalue: the wave- 
function collapses to the eigenvector associated to the 
eigenvalue measured. One can also easily obtain col- 
umn vectors of ISRM by iterating a basis vector. In 
Fig0] we display the Inverse Participation Ratio (IPR) 
for the eigenvectors and column vectors. This quantity 
is given by £ = Yli Wi\ 2 / Yli \^i\ 4 f° r a wavefunction 
l^) = J2i It gives the number of basis states sup- 

porting the wavefunction (£ = 1 for a state localized on 
a single basis state, and £ = N for a state uniformly 
spread over N of them). The results show that where 
intermediate statistics are present, for both eigenvectors 
and column vectors of ISRM one has £ oc N' 1 with 7 < 1 , 
indicating fractal distributions of components. 

In conclusion, we have shown that quantum maps dis- 
playing intermediate statistics can be simulated with a 
remarkable economy of resources on a quantum com- 
puter, especially in a NMR setting. We have also ex- 
plored the link between such intermediate statistics, en- 
tangling power of the quantum evolution and matrix el- 
ements distribution. At last, we have shown that a suit- 



ably randomized map can be used as an efficient gen- 
erator of pseudorandom operators displaying statistical 
properties different from RMT, and in particular pro- 
ducing fractal random vectors. 
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